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Abstract 



o 

^vg ■ It has been speculated earlier that Gamma Ray Bursts are sources of ul- 

trahigh energy cosmic rays. Recently, the search for high energy photons from 
Gamma Ray Bursts by Milagro group has put limits on the isotropic luminosity 
of these transient sources in very high energy photons. The implications of the 
C^^ , results obtained by Milagro to our understanding of the ultrahigh energy cosmic 

^v^ \ ray spectrum from these sources have been discussed in the present work. 
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Gamma Ray Bursts (GRBs) are likely to be excellent sites for particle acceleration 

Ph| and are expected to emit ultrahigh energy cosmic rays (UHECR) [1, 2]. Earlier, 

P ' several groups have investigated the possibility of TeV photon emission from GRBs 

[3, 4] in the burst phase. Mechanisms have been invoked which can give rise to TeV 

^ . energy photons inside GRBs. They are {i) synchrotron self-compton (SSC) process, 

which involves inverse compton scattering of sufficiently high energy electrons on 
'k>( \ the ambient synchrotron photons within the source, the synchrotron photons being 

Vh ' radiated by non-thermal population of electrons inside the GRB, {ii) photo-pion 

production by ultrahigh energy protons within the source and the subsequent decay of 
the neutral pions to photons, {Hi) synchrotron radiation by ultrahigh energy protons 
having energies above 10^°e\^ in the magnetic field within the source. TeV photons 
would interact with the low energy photons inside the source and also with the photons 
of the intergalactic medium. Their flux would be diminished due to their annihilation 
with the low energy photons before reaching the Earth. 

Milagro telescope located in northern New Mexico observes the northern sky for 
lOOGe^ to lOOTe^ 7-ray emission. The Milagro group searched for 40s — 3/i transient 
very high energy emission [5, 6]. They have obtained flux limits and also constrained 
the absolute luminosity of GRBs in very high energy 7-rays. Their limits on very 
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high energy 7-ray luminosity of GRBs is necessarily model dependent. The following 
assumptions have been made in their analysis: (i) the photon spectrum is proportional 
to E~'^ in the energy range lOOGeV^ to 21TeV, (ii) GRBs are isotropically distributed 
in the Universe, [iii) GRBs follow star formation rate in z and (iv) absorption of TeV 
photons by extragalactic background light has been included, following [7] . Assuming 
that all GRBs are emitting a 80sec pulse of very high energy photons, they have found 
that for an average distribution of 4.8GRBs/Gpc^/year over < z < 0.5 the 90% 
confidence upper limit on their isotropic lumisoty is lO^^erg/sec and for an average 
of 0.8GRBs/Gpc^/year over < z < 0.5 the corresponding upper limit on their 
luminosity being 10^'^erg/sec. These limits on luminosity of GRBs also depend on 
the duration of the bursts. Milagro group has used duration of bursts as "80sec" as 
a representative value in their work [6] to illustrate how an upper limit can be set on 
the high energy photon luminosity of GRBs. 

The physical mechanism of very high energy photon production inside GRBs may 
put constraints on the expected UHECR spectrum of protons from these sources. 
The SSC model of high energy photon production does not involve ultrahigh energy 
cosmic ray nuclei and this mechanism has no role in constraining the expected UHECR 
spectrum from GRBs. In the other two models protons are accelerated by internal 
shocks in the magnetic field inside a GRB. The resulting ultrahigh energy proton 
spectrum is expected to be a power law as given below. 
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Ap is the normalisation constant and spectral index a = 2 ± 0.2 [8]. Epb is the break 
energy in the proton spectrum of a GRB. The protons accelerated in an irregular 
magnetic field inside the GRB are emitting synchrotron photons. In the fireball 
model of a GRB, expanding ultrarelativistic shells of matter or wind moving with 
different velocities collide with each other and also with the external medium. If the 
expansion time of the ultrarelativistic shell is longer than the synchrotron energy loss 
time by protons then the protons would loss all their energies by synchrotron photon 
emission, otherwise they would lose only a fraction of their total energies during the 
expansion of the ultrarelativistic wind. The break energy of the proton spectrum 
Epb represents the energy above which the protons would lose all their energies by 
synchrotron radiation within one expansion time scale of the wind. 

Some of the ultrahigh energy protons may come out from the GRBs and give rise 
to the UHECR spectrum observed on Earth [9, 10, 11]. For a = 2, the resulting 
high energy 7-ray spectrum in the photo-pion production model is expected to be 
proportional to Ez^^ in the energy range of our interest. The photo-pion production 



scenario is not an efficient mechanism of very high energy 7-ray production [3], hence 
it has not been included in the present discussion. In the proton-synchrotron model 
of very high energy 7-ray production the 7-ray spectrum is proportional to -E~^'^ and 
E'"^ below and above the photon spectrum break energy respectively for a = 2 [3, 4]. 
Using the proton-synchrotron mechanism of very high energy photon production and 
the limit obtained by Milagro group [6] the contribution of GRBs to the UHECR 
spectrum observed on Earth above energy lO^^eV^ has been explored in the present 
work. 

The Proton- Synchrotron Mechanism of Very High Energy Photon Production: 

This mechanism has been discussed in detail earlier in these two papers [3, 4]. 
GRBs emit 7-rays, X-rays and also photons in optical and radio frequencies. KeV- 
MeV energy photons inside a GRB are likely to be produced by the synchrotron 
radiation of shock accelerated electrons. The energy density of these photons (KeV- 
MeV) in the wind rest frame, f/^, would be almost equal to the energy density of 
the electrons Ue-, if electron synchrotron emission is an efficient mechanism inside the 
ffieball. Initially, the accelerated protons carry nip/me times more energy than the 
accelerated electrons. If we denote the energy density of the accelerated protons by 
Ur, then f/p = — f/g ~ ^^t/^. In this model it is assumed that the magnetic field 

r^ I Trie tyie, ' 

inside the ffieball has an energy density almost equal to the total energy density of 
the ffieball. The total energy density U is nearly equal to Up as most of the energy 
is carried by the shock accelerated protons inside the ffieball. One can relate Ub and 
Up from the above discussion in the following way, 

Ub = ^bUp, (2) 

C,B is the equipartition parameter whose value is of the order of unity. The magnetic 
field energy density can be expressed as, 

f = ^B^U,, (3) 

The energy density of the low energy photons U^ in the wind rest frame can be 
expressed in terms of the luminosity of these photons Ll^ in that frame. The low 
energy photon luminosity of a GRB in the wind rest frame is related to its observed 
luminosity on Earth L°j^ by the factor (1 + z)"^ which takes into account the cor- 
rections in time and energy due to the redshift (z) of the GRB. If the GRB has an 
average Lorentz factor F and the time interval between two successive emission of 
ultrarelativistic shells (variability time scale) as measured by the observer on Earth 
is At"'', then the observed luminosity of that GRB in low energy photons L°A can 



be expressed as a function of the parameters z, F, At"^ and energy density U^. The 
hght curve of a GRB usually shows variation of the intrinsic variability time scale 
in a wide range of values and the intrinsic variability time of a GRB is much less 
than the total duration of the burst. The proton-synchrotron mechanism of high 
energy photon production is not significantly sensitive to the variation in At"** over 
the range 0.1 sec < At°^ < lOsec [3]. The expanding ultrarelativistic shells of matter 
collide with each other and their energies are dissipated. The characteristic distance, 
at which they collide and internal shocks are formed, is called the dissipation radius. 
The dissipation radius as measured at the source in the observer's reference frame has 
been denoted by r^. We would refer to the observer's reference frame at the source as 
the source rest frame where the ultrarelativistic wind moves with Lorentz factor F. 
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where. 



r^ = 2T^cAt°y{l + z) , (5) 

The magnetic field B inside the GRB can be expressed in terms of the GRB pa- 
rameters using equations (3)- (5). The shock accelerated protons inside the fireball 
would emit high energy photons by synchrotron radiation during the expansion of 
the ultrarelativistic shells of matter or wind. The characteristic energy of these high 
energy photons E^ depends on the value of the proton energy Ep and the magnetic 

field B. 

SehBE^ 

where, e is the charge of the proton and the average value of the sine of the pitch 
angle 6 of the proton is n/A. The proton energy and synchrotron radiated photon 
energy can be related as described below using equation (6), 

Ep{GeV) = Ki{E^{GeV)y/^ , (7) 

After expressing B in terms of the GRB parameters we get in the wind rest frame 
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where, F = F300 x 300 and observed luminosity in KeV — MeV photons L°]^ = 
L°^ 51 ^ IQ^^erg/ sec. The energies of any particle in the observer's rest frame and in 
the wind rest frame are related as follows 

E- = ^^E, (9, 



This general relation in equation (9) can be substituted in equation (7) to relate the 
proton energy and synchrotron photon energy in the observer's rest frame on Earth. 

E;\GeV) = K2{E'^\GeV)f'^ , (10) 

where, 

/ts = 2.3 X 10*^^73 300i-^^ (11) 

The fractional energy loss of the protons by synchrotron radiation during one 
expansion time scale (texp) of the ultrarelativistic wind is as follows 

UE,) = J^y (12) 

where, the inverse of the proton synchrotron energy loss time scale is 

' - ' ^"'] (E.), (13) 
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Assuming, Ep >> rripC^ and isotropic distribution of the pitch angles the right side 
of the above equation can be expressed as 
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The expansion time scale is related to the dissipation radius r^ and the Lorentz factor 

r. 

1 r'^ 

The break energy in the ultrahigh energy proton spectrum can be obtained by equat- 
ing the two time scales tps{Ep) and texp- 

Above the break energy Epi, the protons lose all their energies by synchrotron radiation 
during the expansion of the ultrarelativistic wind hence, the fractional energy loss by 
protons can be expressed as 

The expansion time scale texp of the wind and synchrotron energy loss time scale of 
protons tps{Ep) can be expressed in terms of the variability time scale At°^, redshift z 
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and Lorentz factor T of the GRB. Using equation (16) one can obtain the expression 
for the break energy in the proton spectrum in the wind rest frame as given below 
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E„b ~ 1.6 X 10« ' 7 ' — -GeV, 

The break energy in the synchrotron photon spectrum due to the break energy in the 
proton spectrum appears at E^b. 

E,b = ^i"{-^)'GeV, (19) 

■^^ /At°* 



= 1.5 '°° \"; ' GeV, (20) 

In the observer's reference frame on Earth the break energy in the high energy photon 
spectrum is 

8 /At°'' 



e* = 450 ''''" \"' ■' GeV , (21) 

There is an acceleration time scale of the protons ^^tacc" , tacc{Ep) ~ 27n]r L{Ep) / c 
where ^^rL{Ep)" is the Larmor radius, rL{Ep) = Ep/{eB). rj is a factor of the order of 
unity. The shock accelerated protons are losing energy by synchrotron radiation of 
high energy photons. It is possible to obtain the upper cut-off energy in the proton 
spectrum by the condition given below 

The above condition leads us to write the following expression for the maximum 
proton energy in the wind rest frame. 

Y 3/2 /At^y/2 

Ep,ma. = 9.8 X lO^T]-'/^—^ \^ GeV , (23) 

ej'{Lt,y (l + z) 

From equation (9) and (23) the expression for the maximum proton energy to be 
observed on Earth is 
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The maximum photon energy to be observed on Earth can be obtained using 
equation (10) and (24). 



K^:'^ = 1.64 X l0'r^-\-^)GeV , (25) 



We have used 77 = 1 in our numerical calculations. For reasonable ranges of values 
of the GRB parameters, one observes tacc < texp < tps for Ep < Epb < Ep^max, 

T^acc —1 L,ps ^ L,exp 10^' ^'pb ^ ^p _i ^p,rnax) aUQ Lpg <^ t^acc lO^ ^p ^ ^p,max- 

The Ultrahigh Energy Proton Spectrum on Earth: 

The luminosity of a GRB in very high energy photons in the observer's reference 
frame on Earth is 

L^^ = Anir^'T'cUn, , (26) 

, where it has been assumed that the GRB has occured aX z ^ 0. r^ is the character- 
istic dissipation radius measured at the source rest frame. Un-y is the energy density 
of the very high energy photons in the wind rest frame. If the GRB is at z ~ then 
the intergalactic extinction of the high energy photons due to their annihilation with 
infrared background photons can be neglected. If the right side of equation (26) is 
expressed in terms of E^, then 



L'^l = Mr'.fr'c r"''^^'^ E,'^^^^^exp{-nn{E,))dE, , (27) 
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In the present calculation, Eu^^rnin = lOOGeV/F, since Milagro telescope searched for 
photons above observed energy lOOGeV and the upper cut-off energy of the spectrum 
in the wind rest frame is Ewy^max = -^"fmax/r- The internal optical depth Tin{E^) has 
been calculated in [4] for various inferred values of GRB parameters. Above a few 
hundred GeV observed photon energy, the internal optical depth does not change 
with energy. It's average value between EH-y^max and EH-y,min can be used in the 
calculation. We denote it by Cab- The ultrahigh energy proton spectrum in the wind 
rest frame and the high energy photon spectrum expected to be produced by the 
synchrotron radiation of these protons are related as follows 
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Using equation (28) it is possible to write equation (27) in the following form 

< = 47r(r^)^r^c /""^'"^"^ fAEp)Ep^^^^e^,dEp , (29) 
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, where Ep^min and Ep^max can be calculated using equation (7) and equation (23) 
respectively. The ultrahigh energy proton spectrum expected from a GRB has two 
parameters, normalisation constant Ap and the proton spectral index a which are to 
be determined by using the limit on high energy photon luminosity of GRBs obtained 
by Milagro and comparing the ultrahigh energy proton spectrum expected from GRBs 
with the UHECR spectrum observed on Earth by different experiments. Equation 
(29) can be solved to determine Ap in terms of a, L'^^, r^, T, Epi,, Ep^min and Ep 
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where it is assumed that Ep^min < Epb < Ep^max- The number of protons emitted by 
a single GRB per unit observed proton energy on Earth can be written as 

^^'<^"' - tAMrS^/-^ (31) 



dE°>' ' "' dEp 

td is the duration of very high energy photon emission in seconds as observed from 
Earth. Milagro group [6] has used t^ = SOsec for their analysis to show how an 
upper limit can be obtained on high energy photon luminosity of GRBs. Durations 
of GRBs may vary from few hundred milliseconds to tens of seconds. The shortest 
burst detected by BATSE (Burst And Transient Source Experiment) had a duration 
of 5ms and the longest burst so far observed had a duration of dOminutes. One 
may assume an average value for the duration of high energy photon emission from 
all GRBs to obtain an estimate of their contribution to the ultrahigh energy proton 
spectrum to be observed on Earth. The observed proton spectrum on Earth from all 
GRBs is expected to be 
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, where Rgrb is the average distribution of GRBs at 2; ^ 0. dc{E° ) is the distance 
beyond which sources do not contribute to the flux above energy E"''. There is dis- 
cussion on the energy dependent cut-off distance dc{E°^) in [12, 13]. The spectrum 
of particles accelerated by Fermi mechanism in relativistic shocks has been studied in 
detail by several groups [14] . In the ultrarelativistic limit for a pure Kolmogorov mag- 
netic turbulance a spectral index ~ 2.26 ±0.04 is expected [15]. We have obtained the 
expected ultrahigh energy proton spectrum on Earth from GRBs for a = 2.2,2, 1.8. 
In particular, we have investigated about what could be the GRB parameters which 



may give rise to a UHECR proton spectrum on Earth comparable with the observa- 
tional data from AGASA [9] and HiRes [10] experiments above energy lO^^eV. 

Results and Discussions: 

The ultrahigh energy proton spectra expected on Earth from all GRBs have 
been derived for a = 2.2,2,1.8 assuming that only one third of the protons are 
coming out of the GRBs and shown in Fig.l. We have used the following val- 
ues of GRB parameters Lorentz factor F = 300, low energy photon luminosity 
L°j^ = lO^^erg/sec, intrinsic variability time of the source as measured on Earth 
At°^ = 5sec, high energy photon luminosity L°^^ = lO^^erg/sec and distribution 
of GRBs Rgrb = 4:.8GRBs/Gpc^/yr at z ^ 0. We have considered the low en- 
ergy photon (KeV-MeV) spectrum from GRBs with its typical value of break energy 
e^** = 0.5MeV. The spectrum of low energy photons is a power law spectrum in 
energy with spectral indices (3i = 1 and I3h = 2.25 below and above the break en- 
ergy respectively. This low energy photon spectrum has been used to calculate the 
internal optical depth of high energy photons. The internal optical depth has a value 
0.098 at imGeV and its value is 0.196 at lOOOGel/ for V = 300, Lf.^ = lO^^erg/sec, 
At"^ = 5sec, z = 0, ef = 0.5MeV and A = 1, p^ = 2.25. We have used the aver- 
age value of the internal optical depth between observed photon energies lOOGeV and 
lOOOGeV in our calculations. The value of the internal optical depth almost saturates 
to a constant value at higher energies [4]. The particle flux decreases rapidly with 
increasing energies hence, a slight variation in the value of Tin{E^) with increasing 
energies would not affect our results although we are using the exponential of rj„(-E^) 
in our calculations. The break energy, lower and upper cut-off energies in the proton 
spectrum in the wind rest frame are 8 x lO^GeV, 1.67 x lO^GeV and 21.9 x lO^GeV 
respectively. The ultrahigh energy proton spectrum obtained by us for a = 1.8 in 
Fig.l. is comparable with the data recorded by HiRes experiment but below the 
data points of AGASA experiment. Below the peak the spectrum with spectral index 
a = 2.2 is comparable with the data collected by AGASA experiment but it is above 
the data points recorded by HiRes experiment. 

There are several parameters which are involved in this study. It is important 
to dicuss how our results depend on the parameters like the observed variability 
time of GRBs At°^sec, the Lorentz factor F, and the duration of very high energy 
photon emission from GRBs tasec. The internal optical depth {Tin{E^)) of high energy 
photons is inveresly proportional to At°^ and F^ or F'* depending on the photon energy 
[4]. It is directly proportional to the observed luminosity of low energy photons (KeV- 
MeV) Ll\^. 



If At"^ = O.Olsec instead of 5sec, then the value of the internal optical depth 
increases by a factor of 500 and no high energy photons would be able to come out 
from GRBs. An increase in the the value of At"^ to lOsec from 5sec does not give 
rise to any change in the high energy proton or photon spectrum from GRBs. 

Internal optical depth is very sensitive to the changes in the value of the Lorentz 
factor r. If r ~ 1000 instead of 300 and the values of the other parameters remain 
unchanged, then the value of the internal optical depth becomes 123 times smaller 
hence, all high energy photons would be able to come out of GRBs. The break energy 
(Eph) and the upper cut-off energy {Ep^rnax) in the ultrahigh energy proton spectrum 
from GRBs are proportional to F^ and F^/^ respectively. For F = 1000, the break 
energy Epi^ becomes ~ 3.29 x lO^^Gel^ and it is higher than the maximum energy in 
the proton spectrum Ep^ax = 1-33 x 10^'^GeV". In this case the expected ultrahigh 
energy proton flux on Earth from all GRBs would be one order of magnitude higher 
than the observational data of HiRes experiment even after assuming an efficiency of 
proton emission equal to 1/3. If the Lorentz factor is 100 then the internal optical 
depth would be 81 times larger than the corresponding value for F = 300 hence, very 
high energy photons would not be able to come out of the GRBs. The ultrahigh 
energy proton flux on Earth would be very high for F = 100, if we assume that the 
high energy photon luminosity of GRBs is IQ^^erg/ sec as obtained by Milagro group 
for Rgrb = 4:.8GRBs/Gpc^/yr and ta = SOsec. 

If the total energy emitted in very high energy photon emission by each GRB as 
measured on Earth is denoted by E"^^ then, E'^^ = L°^^ x t^. The expected ultrahigh 
energy proton flux on Earth from all GRBs is directly proportional to E'^ . An 
increase in the value of L'^ or td would result in an enhancement of the ultrahigh 
energy proton flux on Earth to be produced from all GRBs. 

Milagro group has assumed the high energy photon spectrum to be proportional 
to -E~^ in deriving the upper limits on high energy photon luminosity of GRBs. We 
have used their limit in the context of the proton synchrotron model of high energy 
photon production. In this model the photon spectrum is proportional to i5~(°'+i)/2 
and ii^~("+2)/2 below and above the break energy respectively where, a ~ 2. [3, 4]. 
The discrepancy in the values of the high energy photon spectral index used in their 
analysis and the present work has to be kept in mind while explaining our results. 
Also, the values of the GRB parameters are not known to us at the time of their very 
high energy photon emission. It may be possible that different GRBs are emitting 
different energies in very high energy photon emission. We do not claim our results to 
be accurate as there are several uncertainities involved in these calculations. We have 
derived what could be the ultrahigh energy proton spectrum on Earth from all GRBs 
if we use the upper limits on their luminosity of very high energy photon emission 
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obtained by using high energy 7-ray telescope. 

Conclusion: 

The ultrahigh energy proton spectra on Earth from all GRBs have been derived 
for a = 2.2,2,1.8 assuming Lorentz factor to be 300, isotropic low energy photon 
luminsoity of GRBs as observed on Earth to be lO^^erg/sec, variability time of GRBs 
to be 5sec as measured on Earth, observed isotropic high energy photon luminosity 
of GRBs to be lO^^erg/sec, duration of high energy photon emission to be 80sec and 
assuming that GRBs are isotropically distributed with a rate of 4.8GRBs/Gpc^/yr 
near z ^ 0. If the efficiency of proton emission from GRBs is assumed to be equal 
to 1/3 then our derived spectrum for a = 1.8 becomes comparable with the data 
collected by HiRes experiment. 
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Figure 1: Proton spectra on Earth from all GRBs have been derived with a = 
2.2,2,1.8 (dotted, dashed, solid lines) and compared with the observational data 
from AGASA [9] and HiRes [10] experiments. The values of the GRB parameters 
used to derive the UHECR proton spectra on Earth from all GRB are as follows: 
r = 300, Lf = 10^^ erg /sec, At"^ = 5sec, L?^^ = 10^^ erg /sec and ta = SOsec. 
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